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We study a network of three populations of coupled phase oscillators with identical frequencies. The popula- 
tions interact nonlocally, in the sense that all oscillators are coupled to one another, but more weakly to those 
in neighboring populations than to those in their own population. Using this system as a model system, we 
discuss for the first time the influence of network topology on the existence of so called chimera states. In 
this context, the network with three populations represents an interesting case because the populations may 
either be connected as a triangle, or as a chain, thereby representing the simplest discrete network of either 
a ring or a line segment of oscillator populations. We introduce a special parameter that allows us to study 
the effect of breaking the triangular network structure, and to vary the network symmetry continuously such 
that it becomes more and more chain-like. By showing that chimera states only exist for a bounded set of 
parameter values we demonstrate that their existence depends strongly on the underlying network structures. 
We conclude that chimeras exist on networks with a chain-like character, which indicates that it might be 
possible to observe chimeras on a continuous line segment of oscillators. 
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Collective synchronization of coupled oscillators 
is a problem of fundamental importance, and oc- 
curs in a wide range of systems including Joseph- 
son junctions arrays, circadian pacemaker cells 
in the brain or the metabolism of yeast cells ^ 
The Kuramoto model has been studied under the 
influence of global (all-to-all) and local (nearest 
neighbor) coupling in great detail. ^"-^^ The in- 
termediate case of nonlocal coupling where the 
coupling strength decays with distance in a net- 
work was first investigated by Kuramoto et al..^^ 
In 2002 they observed a remarkable novel state 
where a population of identical oscillators splits 
into two subpopulations, one being synchronized 
and the other desynchronized, which is called 
chimera state.^^^-^^ Since then, several studies 
have been concerned with its bifurcation behav- 
ior and its emergence under the aspects of het- 
erogeneous oscillator frequencies or delayed cou- 
pling.*^ Chimeras have been observed on a vari- 
ety of network structures such as rings, networks 
with two and three oscillator populations, and 2D 
lattices in the shape of spiral waves. A natural 
question arises: which network structures allow 
for the existence of chimeras? Here we deter- 
mine for the first time limits for the existence of 
chimeras in a simple network of three oscillator 
populations^" as we vary the nature of the nonlo- 
cal coupling amongst the populations. 
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I. INTRODUCTION 

While studying a continuum of identical oscilla- 
tors on a ring with nonlocal coupling, Kuramoto et 
al}^ discovered a remarkable state where the pop- 
ulation of oscillators splits into two subpopulations, 
where one is synchronized and the other is desynchro- 
nized, known as a chimera state. Since then sev- 
eral groups have explored the nonlinear dynamics of 
chimeras4r—i^iS'ii"— 1^124 Their emergence on the ring 
was first analyzed by Abrams and Strogatzp-S: who found 
that chimera states are born through a saddle node bifur- 
cation. The state exists in systems with various network 
structures; for instance, Shima and Kuramoto^ showed 
that chimeras also exist on 2D lattices with free bound- 
aries in the shape of spiral waves. Other authors have 
been studying its appearance on variations of the Ku- 
ramoto model such as systems with delayed couplingj^ii 

While various aspects about the emergence of chimeras 
have been addressed, the question of how the underlying 
network structure affects their existence has not been ad- 
dressed systematically. We make a first step in this direc- 
tion by studying a network of three oscillator populations 
with nonlocal coupling. A network with three nodes may 
either be arranged as a triangle or as a chain, and there- 
fore represents the simplest case of a discrete network 
with ring-like or chain-like character, respectively; the 
network structure will be determined by the nature of 
the coupling, as we explain later. The case of the trian- 
gular network has already been discussed in a previous 
study^° where the author shows that two stable chimera 
states may coexist. In this study, we shall introduce a 
new parameter that allows us to break the rotational in- 
variance inherent to the triangle and to vary the network 



structure continuously such that it becomes more and 
more chain-hke. Thereby we seek to find parameter lim- 
its for the existence of chimera states. 

We consider the case of infinitely many oscillators. 
This case is often considered a valid approximation, and 
using a recently developed method by Ott and Anton- 
seUj ^^i^^ it enables us to reduce this infinitely dimen- 
sional problem to a set of a finite ordinary differential 
equations i^iiS 

The article is structured as follows. In Section II we 
provide the definition of the system under investigation 
and explain how we intend to vary the character of the 
network structure by introduction of a special param- 
eter. We then state the equations resulting from ap- 
plying the Ott-Antonsen method and consider special 
symmetries that allow for chimeras. The analysis of the 
chimera states and their bifurcation scenarios follows in 
Section IIII[ where we analyze how the chimera states 
cease to exist as we vary the network structure using the 
above mentioned parameter. Section fVl summarizes our 
findings. 



II. GOVERNING EQUATIONS 



The governing equations are given by 
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where the phases of the oscillators are defined by 6; i 
denotes the individual oscillators belonging to the popu- 
lation with index cr £ {1, 2, 3}, each of which has N„ os- 
cillators; parameter a changes the character of the phase 
attraction. 

The coupling kernel i^cro-' describes the coupling 
strength between populations a and a' and is assumed to 
decay with increasing separation between the populations 
on the network. Within a population, the oscillators in- 
teract with strength Kcr„> — 1. Neighboring populations 
couple more weakly, with strength 1 — ^orl — c^, as 
displayed in Fig. [TJ We then have 
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In the case where A = we retrieve the case of a glob- 
ally coupled network. Thus A may be thought of how 
'far' we are away from global coupling. The introduc- 
tion of the parameter c in ([2]) generalizes our previous 
workiS: for c = 1, one obtains the case of a triangular 
network with the same rotational symmetry as a contin- 
uum of oscillators on a ring, studied byi^i^iiSiii^ If we let 
c ^ 1 this rotational symmetry is broken; in particular, 
for c > 1, the 'left' and 'right' populations interact less 
strongly with one another than they interact with what 
turns out to be the center population (see Fig. [1] (ni)). 



Hence parameter c controls what one may conceive as 
the distance between the outer left and right population. 
Thus for c > 1, the network attains a chain- like charac- 
ter. To facilitate comparison with earlier work, we focus 
on positive coupling only and impose the constraints 



c- A < 1, 
A<1. 



(3) 
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Our special interest lies in the case where c > 1 because 
of its close relation to chain-like structures. However, we 
relax this constraint during the subsequent analysis in 
favor of a deeper understanding of the bifurcation struc- 
ture. 



A. Reduced equations and symmetry manifolds 

We summarize how we obtain a reduced set of gov- 
erning equations based on the Ott-Antonsen method^-. 
In order to analyze the behavior of ([T]), we consider the 
limit of infinitely large populations, i.e. — > cx). In 
this limit, one may describe the dynamics in terms of 
the oscillator density distribution f'^{9,t), which evolves 
according to the continuity equation 
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The velocity of the oscillators in this limit is given by 



3 

E 



sm 



{9' -e-a)r'{e',t)de\ 



and a complex order parameter is introduced. 
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defining an average over all oscillator phases; it is there- 
fore referred to as weighted or global order parameter. In 
this formulation of the problem, we are able reduce the 
problem to a finite set of equations using a method re- 
cently developed by Ott and Antonsen.^* Following Ott 
and Antonsen, we restrict our attention to a special class 
of density functions in the form of 
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where c.c. is the complex conjugate of the expression un- 
der the sum. The application of this ansatz ultimately 
leads to a set of ordinary differential equations. How this 
reduction method is performed on the problem of oscilla- 
tor populations and leads to a finite set of equations has 
been reported in detail elsewherej ^'^^i^^ We skip the full 
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FIG. 1. Resulting network structures by varying parameter c. The gray disks symbolize populations, inhabited by individual 
oscillators symbolized by black dots. Their bidirectional coupling is represented by black lines. Each population has a self- 
coupling of unit strength 1. The population in what becomes the center for c 7^ 1 is coupled to the neighboring populations 
with strength 1 — A; the populations to the left and right are coupled with strength 1 — cA. The case of a triangular network 
is obtained for c = 1; the character of the network has chain-like character for c > 1. 



derivation and state the resulting reduced equations: 

Pa = ^ Kca'Pa' sin {(jic,' - (j>a + 13), (9) 



where 
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/3 = 7r/2 - a, 
aa = Pae^"''", 
a e {1,2,3}. 
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The variable a* takes account for the local order for each 
population and found to be related to the gfobal order 
parameter via 
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The dynamics of the oscillators is described by two vari- 
ables per population: p„ describes their degree of syn- 
chronization and 4>a represents the average phase of the 
oscillators in each population. 

Before analyzing these equations, we restrict our at- 
tention to specific symmetry manifolds. Perfectly syn- 
chronized populations have p^- = 1 and desynchronized 
populations have po- < 1; the latter consists of oscillators 
that drift relative to one another and to the synchro- 
nized populations. Let S and D denote synchronized and 
desynchronized populations, respectively. In the trian- 
gular network, due to its rotational invariance, we could 
only distinguish two chimera states, namely SDS (sync- 
drift-sync) and DSD (drift-sync-drift). The situation is 
different here, and other chimeras are possible in the case 
of c 1: SDD, SSD and their symmetry-equivalent re- 
flections across the center population. These states are 
however not the focus of this study and are excluded in 
our analysis. 



The SDS state is defined as pi = p3 = 1 and p = 
P2 < I, whereas the DSD state has p = pi ^ p^ < I and 
P2 = 1. Because our coupling kernel ^ is symmetric. 
Pi — P3 also implies 0i — (pa; hence populations 1 and 
3 are phase- locked (this holds true also if c 7^ 1.). The 
phase difference of the angular order parameter between 
the synchronized and desynchronized states is defined by 

Ip ^ (pi - (1)2 = (1)3 - Cp2- (15) 

Applying these symmetry assumptions to ^ and substi- 
tuting the coupling kernel defined in ([2]), we obtain the 
equations describing the 5*1? 5' states 

p = — ^[2(1-^) sin(^ + /3)+psin/3], 
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Fixed points of ((T6| and (jl7l) correspond to phase-locked 
solutions of the original system. The reduced equations 
© (after transformation into the co-rotating frame, i.e. 
(p^ ^ (p^ -\- ujt) , and (jl6ll7p share the property of being 
invariant under the following time-reversibility transfor- 
mations: 



(/3,i,V)^ (-/3,-i,-V') 
(/3,t)-^ (/3 + 7r,-t). 



(18) 
(19) 



The bifurcation structures discussed in the next section 
repeat themselves accordingly in the (/3, A)-plane. 

In conclusion, we have reduced the governing equa- 
tions ID) to a low dimensional system for the local order 
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parameters, and we have cast our problem into a two 
dimensional system represented by Eqs. (jl6p and (TTTl) . 
This enables us now to study the problem in the phase 
plane. 

III. BIFURCATION BEHAVIOR NEAR THE 
TRIANGULAR STRUCTURE 

We briefly review the findings for the triangular cou- 
pling (c = 1, before we study the symmetry breaking 
case (c 7^ 1) in more detail. The associated bifurcation 
diagrams are shown in Fig. [31 where we compare the 
triangular case (left) with the symmetry breaking case 
(right). We first consider the case of SDS symmetry. 
We keep the parameter /3 constant while increasing the 
value of A step by step: close to global coupling, i.e. for 
small values of A, we only observe the always present in- 
phase SSS solution. As we increase A, a stable chimera 
is born in a saddle-node bifurcation. One step further, 
the chimera loses its stability through a Hopf bifurca- 
tion and a limit cycle is born, which corresponds to a 
so-called breathing chimera with an oscillating order pa- 
rameter. Increasing A further, we observe how the limit 
cycle grows and eventually collides with the saddle and is 
destroyed in a homoclinic bifurcation. The saddle-node 
curve, Hopf and homoclinic bifurcation curves all inter- 
sect in the Bogdanov-Takens (BT) point of codimension 
two. The DSD symmetry exhibits a similar bifurcation 
structure: for small A, the scenario is identical to the one 
seen for the SDS symmetry; however for larger A, the 
scenario is repeated in reversed order (with increasing 
values of A), as shown in Fig. [5] 

We consider now the symmetry breaking case relatively 
close to the triangular symmetry with c — 1.1. The bi- 
furcation diagrams are obtained by inspection of phase 
portraits (large dots in the figure) . Saddle node and Hopf 
curves may be calculated analytically by simultaneously 
solving the fixed point equations (fT6)) and (flT)) with ei- 
ther the saddle node condition, 

det (J) = (20) 

or the Hopf condition, 

tr( J) = and det (J) > 0, (21) 

where J is the Jacobian of (|16l) and (IT71) . The resulting 
equations are solved using a series approach with /3 as a 
bifurcation parameter, similar to the method described 
in detail in earlier work,-- however by assigning a desired 
value to c before computing the series coefficients. The 
perturbation series for A — A(jS) yield no further insight 
and we omit them here for brevity. We observe a similar 
bifurcation scenario as before in the triangular case, but 
there is an important qualitative change: the parameter 
region allowing for SDS chimeras (above) has shrunk in 
parameter space, whereas the attractor region of DSD 
chimeras (below) has grown. Interestingly, it looks as if 
the homoclinic curves for the DSD chimera now coincide 



for a wide range of /? values. Unfortunately we could not 
confirm this behavior analytically, as we were unable to 
determine the Melnikov integrals leading to homoclinic 
conditions, but the phenomenon is not essential for the 
matter of this study. 

Moreover, the Bogdanov-Takens (BT) points have 
moved to the left and right, respectively, according to 
the shrinking or growing of the regions of existence for 
chimeras. We can check this by determining the locus of 
the BT points for c = 1.1. We numerically solve the fixed 
point equations implied by Eqs. (jl6p and (1171) simultane- 
ously with the saddle node ((20)) and Hopf condition (j2ip . 
For c = 1.1 we find the following critical points for the 
SDS symmetry, 

(/3,A)si?s« (0.0902,0.5361), (22) 

and for the DSD symmetry, 

(/?, A)dsd,i « (0.2467, 0.6466) (23) 
(/3, A)osD,2 « (0.2244, 0.8132), (24) 

which are denoted by arrows in Fig. [2l Comparison of 
these numbers with our previous results^*', for c = 1 con- 
firms our observation that the BT points have shifted to 
the left and right, respectively. We proceed with finding 
analytical expressions describing how these points move 
in parameter space and put limits on the existence of 
chimeras. 



IV. LIMITS OF EXISTENCE FOR 
CHIMERAS 

We study now more generally what happens to the 
chimeras as we vary the values of c. Bifurcation diagrams 
for SDS and DSD chimeras were obtained by numerical 
continuation for a range of c values (Fig. [3]). The behavior 
described above that the attracting region for the SDS 
chimeras shrinks and the region for DSD chimeras grows 
with increasing values of c, is now confirmed for a larger 
range of parameter values; conversely, as we consider the 
case c < 1, the opposite holds true. Regions of existence 
for chimeras disappear completely at the point where the 
BT point touches the /3-axis: this is where the saddle 
node, Hopf and homoclinic curves collapse in a single 
point. 

We calculate the critical c values where this occurs. 
The /3-axis represents a singular limit and simply sub- 
stituting /? = into the equations defining the BT 
points (saddle node ((20)) and Hopf condition ((2T))). leads 
nowhere; however we can calculate Ccrit by using a per- 
turbative approach. We use /3 as a perturbation param- 
eter in the limit /3 — 0. Inspecting the phase portrait of 
(IT6l) we find that the SDS chimera is located near -0 = 0; 
the correct perturbation ansatz is 

= 0i/3-f 0(/32). (25) 
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FIG. 2. Effect of breaking the rotational symmetry on the bifurcation diagram for the SDS and DSD symmetries. The 
triangular case-^ with c = 1 is shown in the left column in (a) and (c) for comparison with the case of broken symmetry 
(c = 1.1) in the right column in (b) and (d). The displayed curves are: the saddle-node curve (solid red), the Hopf curve 
(dashed blue), and the homoclinic curve (dotted black). Dots mark the bifurcation points obtained by inspection of the phase 
plane. The homoclinic curve is an interpolation based on these points, whereas the solid curves were obtained analytically. 



Solving the resulting equations at first order of /?, we find 
that the critical point for the SDS symmetry is located 
at 



^ (l5 - \/33 + 2\/6\/33- 30j w 1.12356, 
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0.770111. 



(26) 



By the same token, for the DSD chimera associated with 
small values of A (lower DSD chimera), we have 

c^„t = (a77A - 350%/57 + (473%/7 - 5SVS99j 



X 



\Jil + V57j /3976 « 0.677869, (27) 



l + \/l + y57 



399 - IblVf) /1988 « 0.735569,(28) 



Pcrrt = 7 ( 1 



0.552586. 



(29) 



Observing Fig.|31 we notice that the saddle node curves 
of the upper DSD states coincide for all values of c; hence 
the associated set of BT points must also be located on 
this curve. Following the motion of the BT point as c de- 
creases, we suspect that c becomes critical exactly when 
the BT point has reached A = 1. (This looks like a limit 
of a higher order singularity than previously; indeed it is 
much harder to see in phase portraits at which angle ip 
the BT point detaches from the /3-axis.) This time the 
correct perturbation ansatz turns out to be 



^ = --fVi/3 + 0(/?2). 



(30) 



Solving again to first order in /3 results in 
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FIG. 3. Bifurcation diagram for the SDS chimera (above) and 
the two DSD chimeras (below) in the (/3, yl)-plane for a range 
of c values. Saddle node (solid red) and Hopf curves (dashed 
blue) are shown. The curves related to the SDS chimera 
collapse onto the /3-axis at c « 1.123; conversely, the two 
DSD curves collapse on the axis at c « 0.6778 and c = 2/3, 
respectively. It is seen that the Bogdanov-Takens point (black 
dots) of the upper DSD chimera follows the associated saddle 
node curve as we vary c. 



Pent = -y^ • (31) 

We conclude that SDS chimeras exist for parameter val- 
ues c < 1.12356, and DSD chimera for c > 0.677869 
and for c > 2/3, respectively (we shall discuss below that 
these limits for the existence of stable chimeras hold if 
/3>0). ^ 

Another interesting phenomenon is seen in Fig. |31 the 
saddle node curves of the SDS chimera pass through 
the origin (/?, A) = (0,0) only if c < 1, but detach from 
the origin as soon as c > 1. The analogous behavior is 
seen for the DSD chimera, however for the converse case 
where c < 1. Thus having passed c = 1 and approach- 
ing Ccrit implies the shrinking of the chimera attractor 
regions not only in the direction of the /3-axis, but also 
along the direction of the A-axis. To shed more light on 



this behavior, we calculate the locus of the saddle node 
transition A = A{c) for /3 — 0. These curves are this 
time simply obtained by letting /3 = while solving the 
saddle node condition together with the fixed point con- 
ditions. For SDS, we have: 

A = ^ (l2c2 - 9 ± V3v/(3 - 2c)3(l + 2c)) , (32) 

A = 0, (33) 

and for DSD, we obtain the equations 

27(yl -2){A- l)2yl = Ac [-54 -I- 18^(6 + c) 
+A^c {18 + c^) 
- 2A\27 + c{18 + c))] , (34) 
A=l. (35) 

The resulting curves are shown in red in Figs. 2] (a) and 
(b). Recall that we are only interested in regions repre- 
senting positive coupling, and regions where the coupling 
decays with increasing distance on the network graph. 
The second and forth quadrants in the (c, A)-plane rep- 
resent regions where the coupling matrix has a non- 
decaying character because of its entry with (1 — cA) > 1 
in and has therefore no chain-like character. The en- 
tire lower half-plane has the same problem, but now due 
to the coupling kernel entry {I — A) > 1. In the quadrant 
that is left, we have to respect the constraints A < 1 and 
1 < c ■ A; they are indicated by gray dotted curves. All 
these regions that are not of interest to us are shaded in 
gray. The curves defined by ((32)) -(|35 p reach into these 
forbidden areas and are related to interesting dynamics; 
its nature is beyond the scope of this study and would 
require a more in-depth analysis. 

Special points of interest are labeled as follows: A is 
where the saddle node curve detaches from the /3-axis, 
as shown in Fig. [3J the various 's represent the criti- 
cal points where the BT points, together with their sta- 
ble and breathing chimera regimes, collide with the (3- 
axis, for SDS and DSD symmetries; C denotes the point 
where the 'upper' DSD chimera boundary psp intersects 
with the 1 = c ■ A boundary; and finally, D denotes the 
analogous intersection with the A = I boundary. 

Note that the curves described by Eqs. do not 

delineate complete regions of stability: they represent the 
boundary of lowest possible A values where a saddle-node 
bifurcation may occur for at least one /? (this point always 
lies on the (3 = axis due to the observed monotonicity 
of the saddle- node curves). In other terms, chimeras that 
have been created through a saddle-node transition while 
increasing A may lose their stability through a Hopf and 
get destroyed in a homoclinic bifurcation for larger values 
of A, depending on the value of /3. Thus the curves should 
be understood as lower boundaries (or upper for A = 1, 
respectively, for the second DSD state) where chimeras 
still exist. Accordingly in Fig. 21 regions where chimeras 
cannot exist at all are labeled with No Chimera, and 
regions where stable chimeras exist (for an appropriate 
choice of /3) with Chimeras. 
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FIG. 4. Boundaries for the occurrence of saddle node tran- 
sitions in the (c, yl)-plane (i.e. saddle node curves at /3 = 0), 
shown for the SDS symmetry (a) and the two DSD symme- 
tries in (b) (A = 1 is the boundary for the second DSD 
state seen in the upper part in Fig. [3] b)). The regions 
shaded in gray either have negative coupling or have cou- 
pling without chain-like character, as explained in the text. 
Throughout the regions labeled with No Chimera, we find no 
chimeras. The black dots indicate points of special interest: 
A: A = (1,0) is where the saddle node curves detach from 
the origin (/3, A)-plane. B: are the points for which the BT 
points collide with the /3-axis leading to the annihilation of 
the chimera state; B « (1.123,0.543), Bi ^ (0.677,0.735) 
and B2 ~ (2/3, 1). C: intersection with boundaries of pos- 
itive coupling; Csds = ((27 -I- ^)/26, 26(27 + V27)) and 
Cdsd = (1, 1). D: intersection of the saddle-node boundary 
with 1 = A; D = (3/2, 1). The regions where chimera exist 
for /3 < 0, which are stable within the symmetry manifolds 
SDS and DSD, are hatched green. 



Finally, we have a brief look at states with f3 < which 
so far have been neglected in previous studies. The fol- 
lowing applies for both SDS and DSD chimeras, but for 
simplicity let us consider only the SDS symmetry. At 



point B, the BT point collides with the /3-axis, and we 
expect not to see any chimeras beyond this point. But 
for c > Ccriti after crossing the curve defined by BD, one 
observes the creation of a saddle and a source with p < 1, 
corresponding to an unstable chimera. However, virtue 
to the symmetry in p^ . {P,^,t) — > (— /3, — -0, — t), a re- 
lated stable node exists for negative /3. Hence if we ex- 
tend our study to include /3 < 0, stable SDS chimeras are 
also found for c > Ccrit^ as well as the two DSD chimera 
states beyond their corresponding Ccrit-values. The green 
hatched areas denote their regions of existence for /? < 0. 
Still, they leave our domain of interest in Point C and 
completely cease to exist in D. We note that no Hopf or 
homoclinic bifurcations occur in these regions: they ap- 
pear to be completely annihilated in point B (where the 
BT points collide with the /3-axis). Moreover, accord- 
ing to the symmetry defined by there should also 
be instances of stable chimeras for c < Ccrit and /3 < 0: 
the chimeras that are stable or breathers for /3 > are 
present with inverted stability, i.e. the source inside the 
limit cycle becomes a stable chimera for /3 < . 

It is however conceivable that these states with /3 < 
are unstable towards perturbations in the transverse di- 
rections of the symmetry manifolds SDS and DSD (just 
as it seems to be the case for the DSD in the upper part 
of the {A, /3)-planej^) their occurrence is however beyond 
the scope of our interest and we leave this question for 
others to solve. 



V. DISCUSSION 

We have investigated the problem of three nonlocally 
coupled oscillator populations, thereby generalizing work 
on networks with triangular symmetryJ^ By introducing 
the structural parameter c, we were able to study the 
behavior of oscillator populations as we change the qual- 
ity of the nonlocal coupling kernel ([2]). These changes 
differ from the quality changes induced by parameter A, 
controlling how close we are to local or global coupling. 
Conversely, parameter c controls the distance between 
the 'left' and 'right' populations, and we may distinguish 
the different qualities it tunes into as visualized in Fig. [1] 
(i) for very small c, the left and right populations are 
very close, and they see the center as a 'satellite' popu- 
lation. In this limit, the presence of the SDS chimera 
is dominating. One might argue that it corresponds - 
in our terminology - to the SD chimera observed in two 
oscillator populations by Abrams et al.^; (ii) for c = 1, 
we retrieve the triangular case with rotational invariance; 
(iii) for c > 1 the network acts like a chain (discussed be- 
low) and (iv) for very large c, the outer populations will 
almost only sense each other's motions indirectly via their 
coupling to the center population, as long as 1 > c • A, 
i.e. the coupling stays positive. 

We determined lower and upper limits (OBI) for c, (P7)) 
and ([3T|). where chimeras cease to exist. We found that 
these limits are valid for /3 > 0, but that stable chimera 
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may exist (within other hmits) beyond the critical c val- 
ues for /3 < 0. Furthermore, we were able to find limit- 
ing curves A = A(c) defining coupling parameter regions 
(c, A) for which chimeras cannot exist at all. 

We found that parameter c needs to stay relatively 
close above 1 to observe a chain-like SDS chimera state. 
Regardless this closeness to the triangular symmetry c = 
1, any case with c 7^ 1 may be considered a valid instance 
of the discretization of a chain. For example, consider 
the case of a line-segment subject to a coupling kernel 
with exponential decay, e"*^'^'; then we may choose the 
characteristic length scale of the kernel to match a 
desired value of c. It is therefore likely that chimeras may 
also exist for a continuum of oscillators on a line segment. 

It is worth noting how the effects of the coupling pa- 
rameters A and c are related to coupling kernels studied 
in other chimera systems. Abrams et al.^ discussed the 
emergence of chimeras on a continuum of oscillators on 
a ring and thereby used a kernel G{x) = 1-1- A cos (a;). 
Clearly, the effect of A is the same as in our study by 
modulating the amplitude of the coupling. In an ear- 
lier study (also on a ring), Kuramoto and Battogtokh^^ 
used an exponential kernel G{x) ^ exp{—K\x\) where n 
tunes the width of the nonlocal coupling. If space is 
discretized, the effect of varying this width is similar to 
controlling the 'distance' between the resulting oscilla- 
tors populations, just like our parameter c does it here. 
In this sense, our choice for the coupling unites both qual- 
ities. In all systems using kernels with a modulation pa- 
rameter A, it is found that chimeras emerge in the limit 
where /3 = 0, where the system is nearly conservative. 
On the other hand, in a recent study^^ we found that 
spiral wave chimeras in the plane - coupled via a Gaus- 
sian G{x) = exp (— |a:p) - emerge from the different limit 
/? = 7r/2 (i.e. near the gradient system limit). The pre- 
cise origin for this adverse behavior or how the emergence 
of chimera states in these quite distinct limits is related 
to the different types of coupling remains to be solved. 

The issue of which classes of networks allow for chimera 
states to exist in general is an open problem. A related 
matter is the characterization of the basins of attraction 
leading to chimera and globally synchronized states, a 
question raised just recently.-12 For future research, we 
suggest to do more studies along these lines. 
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